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1.0 Executive Summary 


A goal of the National Aeronautics and Space Administration (NASA) Earth Science 
Enterprise (ESE) Earth Science Applications Program is to infuse NASA remote sensing 
data sources into existing partner agency decision support tools in order to enhance the 
performance of these tools. Through IDEA (Infusing satellite Data into Environmental 
Applications) NASA, in partnership with the U.S. Environmental Protection Agency 
(EPA), has performed this data enhancement on a project, in which NASA data was 
utilized to improve particle pollution forecasts. Researchers from NASA Langley 
Research Center and the EPA used data from the Moderate Resolution Imaging 
Spectroradiometer (MODIS) sensor and combined it with EPA ground network data to 
create a NASA data enhanced Forecast Tool. This tool is used to assist forecasters with 
providing the forecasts of particle pollution, or particulate matter less than 2.5 microns in 
diameter (PM 2 . 5 ), for the EPA’s Air Quality Index (AQI). The goal of this project is to 
use existing data sets and models developed for tropospheric chemistry research to aid 
the EPA and state and local agencies in making decisions concerning air quality to 
protect public health. 
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2.0 Introduction 


The Earth System responds to both naturally occurring and human-induced change. The 
National Aeronautics and Space Administration (NASA) Earth Science Enterprise (ESE) 
seeks to understand the response of the Earth System via long-term observations from 
ground networks, sub-orbital platforms, and space-based assets. The role of the Earth 
Science Applications (ESA) Program within the ESE is to incorporate these observations 
into decision support tools employed by partners and to assess the performance of these 
measurements in decision support tools. The approach is to enable the incorporation of 
Earth Science mission outputs (i.e., models and remote sensing data products) to serve as 
inputs to decision support systems. Ultimately, the desired outcome is an enhanced 
decision support tool that results in significant socio-economic benefits. 

One application into which NASA observations have already been incorporated is Air 
Quality Management. During the fall of 2003, NASA, through the Infusing satellite Data 
into Environmental Applications (IDEA) project, provided a prototype, near real-time 
data-fusion product to the Environmental Protection Agency (EPA) with the goal of 
improving the accuracy of EPA’ s next-day Air Quality Index (AQI) forecasts, (Kittaka, 
2004; Szykman, 2004; and Al-Saadi, 2004). 

2.1 Criteria pollutant and scientific rationale 

Under the Clean Air Act of 1990 (http://www.epa.gov/oar/caa/contents.html), the EPA is 
required to set standards for concentrations of air quality pollutants, ensure these 
standards are met through monitoring, and establish a consistent means of reporting air 
quality to the public, which, currently, is the Air Quality Index (AQI). The EPA is 
currently setting air quality standards relating to the concentrations levels of six main air 
pollutants: ozone, particulate matter, carbon monoxide, nitrogen dioxide, sulfur dioxide 
and lead. The EPA is also responsible for forecasting concentration levels of two of these 
pollutants: ozone and particulate matter (PM). These forecasts are used to alert the 
public about potentially harmful ozone and PM levels. Particle pollution or particulate 
matter is the general term used to describe a mixture of solid particles and liquid droplets 
in the air. Particles are classified as coarse (i.e., dirt or dust) or fine (i.e., the by-products 
of fuel combustion). PM can be emitted directly into the atmosphere (i.e., wind blown 
dust or dirt from unpaved roads) or formed in the atmosphere through chemical reactions 
(i.e., sulfates and nitrates formed from emissions from power plants and vehicle exhaust.) 
Many voluntary programs exist between the EPA, state agencies and industry, which 
enable “Action Days” once pollution levels exceed a certain threshold. In addition to 
notifying the public to limit exposure, voluntary emission reductions are suggested for 
industry and for private citizens. 

In October of 2003, the EPA began providing AQI forecasts for particulate matter less 
than 2.5 pm in diameter, or PM 2 . 5 . Particulates in this size range are called respirable 
aerosols and are easily entrapped by the lungs. Pollutants and diseases carried by 
respirable aerosols are a significant health threat. According to the World Research 
Institute, an environmental research and policy organization, “the health effects of 
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particulates are strongly linked to particle size. Small particles, such as those from fossil 
fuel combustion, are likely to be the most dangerous, because they can be inhaled deeply 
into the lungs, settling in areas where the body’s natural clearance mechanisms can’t 
remove them” (WRI, 1999). The EPA notes that the chief causes for concern are that 
increases in PM levels are linked to: 

1 . Increased hospital admissions and emergency room visits for people with heart 
and lung disease 

2. Increased absences from work and school 

3. Reduced visibility due to haze 

4. Altered nutrient balance in the soil and in bodies of water where PM settles 

5. Stained and/or eroded buildings, historical monuments, etc. which are costly to 
repair 


Index 

Values 

Category 

Cautionary Statements 

pm, 5 

(ug/m 3 ) 

PMi 0 

(ug/m 3 ) 

0-50 

Good 

None 

0-15.4 

0-54 

51-100 

Moderate 

Unusually sensitive people 
should consider reducing 
prolonged or heavy exertion 

15.5-40.4 

55-154 

101-150 

Unhealthy for 
Sensitive Groups 

Sensitive groups should reduce 
prolonged or heavy exertion 

40.5-65.4 

155-254 

151-200 

Unhealthy 

Sensitive groups should avoid 
prolonged or heavy exertion; 
everyone else should reduce 
prolonged or heavy exertion 

65.5-150.4 

255-354 

201-300 

Very Unhealthy 

Sensitive groups should avoid 
all physical activity outdoors; 
everyone else should avoid 
prolonged or heavy exertion 

150.5-250.4 

355-424 


Table 1 . EPA Air Quality Index for Particle Pollution. 

In early 2003, the EPA received a Congressional mandate to revise the standards that 
govern PM2.5. The standards being revised included peak concentration values per hour, 
peak concentration values over a 24-hour period for several different land cover types, 
and compliance penalty thresholds relating to litigation and public health. Technical 
definitions of legal requirements were revised and preparations for monitoring 
intercontinental transport were made. The timing of this mandate, along with the 
development by the EPA of standard methods for forecasting PM2.5 and the availability 
of a mature satellite aerosol product at relatively high spatial resolution, made PM2.5 
forecasting an ideal candidate for incorporation of NASA satellite data. 
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The EPA has a ground network of monitoring stations around the country that are 
currently being used to monitor concentration levels of PM 2.5 and as input to forecasts of 
PM 2 . 5 . However, the EPA recognizes that the utility of this network may be improved 
with the addition of satellite imagery, which can provide information about the air in 
regions not covered by these monitoring stations. The addition of NASA data could also 
assist in identifying areas that are generating particle pollution and areas that are 
receiving pollution due to transport between regions. Aerosol products from existing 
NASA satellite systems were identified as potential data sets that could add value to the 
forecast. 

2.2 Ground based measurement characteristics 

The EPA ground network consists of in-situ PM 2.5 monitoring stations located throughout 
the country, with higher concentrations of monitors in more densely populated regions. 
These monitors are operated by the State and Local Air Monitoring Stations (SLAMS) 
and National Ambient Monitoring Stations (NAMS) networks. In addition to Federal 
Reference Monitors (FRM), which acquire measurements of PM 2.5 over 24-hour sample 
periods, the ground network consists of several hundred continuous PM 2.5 Federal 
Equivalent Monitors (FEM) that report PM 2.5 concentration data hourly to the EPA’s 
AIRNow Data Management Center. The continuous FRMs utilize several methods that 
measure the different properties of suspended particles, including mass and mass 
equivalent, and visible light scattering. The Tapered Element Oscillating Microbalance 
(TEOM) instrument measures particle mass as determined by its inertia, with a detection 
limit of ~5pg/m 3 for a 5 minute average. The Beta Attenuation Monitor (BAM) measures 
particle mass by its electron attenuation properties, with a detection limit of ~5pg/m 3 for 
a one hour average. The nephelometer measures light scattered from particles and gases 
and provides a direct estimate of the aerosol light-scattering coefficient with a detection 
limit of ~Mm ' 1 for a ten minute average. Further details can be obtained from Watson et 
al, (1998). 


2.3 Satellite based measurement characteristics 
2.3.1 MODIS 

The Moderate Resolution Imaging Spectroradiometer, MODIS, aboard the Terra satellite 
has a 10:30am equatorial overpass time along a sun-synchronous near-polar orbit. 

MODIS has a viewing swath width of 2,330 km and provides global coverage every one 
to two days. For additional information regarding MODIS, refer to 
http : //modis .gsfc .nasa. gov. 

The MODIS data products utilized in this study are the aerosol optical depth (AOD) and 
the cloud optical thickness (COT). Aerosol optical depth retrieved at ,55pm, which is 
one of the parameters included in the aerosol data product, MOD04 L2, has a resolution 
of 10x10 km 2 at nadir (Kaufman et al, 1998). This aerosol product provides a measure of 
extinction (how much light is unable to pass through a column of atmosphere as a result 
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of aerosols, or particles, in the air) and therefore can be used to estimate the amount of 
aerosols in the atmosphere. Cloud optical thickness included in the cloud product, 
MOD06L2, has a resolution of lxl km 2 at nadir. The COT provides cloud locations as 
well as cloud radiative properties. In this study, the COT was degraded to 5x5 km 2 for 
visualization. 

2.3.2 GOES WF_ABBA 

The Geostationary Operational Environmental Satellites (GOES) satellites are part of the 
NOAA operational weather satellite system. They are in geosynchronous orbits allowing 
them to maintain fixed positions relative to the Earth. Additional GOES information can 
be obtained from 

http://orbit-net.nesdis.noaa.gov/arad/fpdt/goescat_v4/html/GOES_I_l_overview.html. 

Wildfire locations are determined from the Wildfire Automated Biomass Burning 
Algorithm (WFABBA). The WFABBA data is courtesy of NOAA/NESDIS (National 
Environmental Satellite, Data and Information Services) and the University of Wisconsin 
Cooperative Institute for Meteorological Satellite Studies (CIMSS). The GOES 
WF_ABBA is an automated algorithm which uses visible and infrared wavelengths to 
locate fire pixels and characterize sub-pixel burning. General information for 
WF ABBA is located at http://cimss.ssec.wisc.edu/goes/bum/abba.html. Detailed 
information on the GOES fire algorithm and recent improvements can be found in Prins 
and Menzel, (1994) and Prins et al., (1998; 2001a; 2001b). 

2.3.3 Eta/EDAS 

The Eta Regional Forecast model and associated Eta Data Assimilation System (ED AS) 
is part of the operational forecasting system within the National Weather Service's 
National Centers for Environmental Prediction (NCEP). Details of the NCEP forecasting 
and data assimilation system can be found at http://wwwt.emc.ncep.noaa.gov/modelinfo/. 
Eta forecasted winds are used to provide trajectory-based 48-hour forecasts of the 
movement of regions of high MODIS AOD to the EPA PM 2.5 AQI forecasters. The Eta 
forecast winds are provided in GRIB format and are obtained at ftp://tgftp.nws.noaa.gov. 
The NOAA Air Resource Laboratory (ARL) provides an archive of the NCEP ED AS 
analyses. For this study, the u- and v-components of the wind at the 850 mb pressure 
level were used to provide a meteorological context for the MODIS and PM 2.5 
observations in the IDEA data fusion product (see Figure 2). Additional ARL EDAS 
archival information is available at http://www.arl.noaa.gov/ss/transport/edas.html. 

2.4 Objective of comparison 

The objective of the comparison is to determine the appropriateness and benefit of 
additional forecasting tools derived from combining MODIS fine aerosol data and 
modeling code developed for tropospheric chemistry research at NASA Langley. The 
forecasting tools are provided to Air Quality Index forecasters through the EPA to 
potentially improve the forecasting of fine aerosol pollution outbreaks. 
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3.0 Site-by-Site Satellite/In-Situ Comparison 

3.1 Background on time period 

The time period chosen for the comparison is September 2003. Evaluation is needed for 
a PM 2.5 pollution event and September marks the peak in the forest fire burning season. 
Figure 1 shows a map of fire locations during September 2003. The pink and purple 
diamonds indicate a 24-hour accumulation of fire locations as detected by the GOES 12 
Wildfire Automated Biomass Burning Algorithm (WFABBA). Pink indicates the 
presence of a fire whereas purple indicates high probability of a fire location. For 
additional information about GOES WF-ABBA refer to 
http://cimss.ssec.wisc.edu/goes/bum/abba.html. 


20030901-20030930 WF ABBA fire counts 



Figure 1. GOES 12 Wildfire Automated Biomass Burning Algorithm (WF ABBA) fire 
locations for September 2003. 

During September there were frequent small fires associated with agricultural burning in 
the lower Mississippi River valley and Alabama. In early September 2003 forest fires in 
the Northwestern United States and British Columbia produced emissions that led to a 
large enhancement in tropospheric aerosol loading. The initial aerosol loading from these 
Northwestern wild fires in the troposphere was captured by MODIS. This enabled the 
IDEA team to produce forecasting tools during a pollution outbreak and evaluate their 
effectiveness and impact on air quality forecasting in the Midwest and Eastern United 
States. 
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Figure 2. Satellite/in-situ fusion map for 1 September 2003. 

Figure 2 is an example of one forecast tool produced by the IDEA team. It depicts a 
composite integrating data from several sources on September 1, 2003 (indicated at top of 
figure). Daily maps of this fusion product from September 1 through September 30, 2003 
are provided in Appendix A. Aerosol optical depth (AOD) from the MODIS instrument 
aboard the Terra satellite is shown on the rainbow color scale with values greater than 1.0 
shown in dark red. MODIS cloud optical thickness (COT) is shown on the gray scale. 
MODIS AOD is not derived for cloud filled pixels, therefore the cloud fields derived 
from MODIS are plotted to help define synoptic weather features important to the 
movement of aerosols. The location of each continuous PM2.5 monitor from the SLAMS 
and NAMS is represented by a colored circle. The color of the circle indicates an hourly 
averaged PM2.5 concentration at 20Z with a color scale shown in the vertical color bar. 

The US EPA Air Quality Index (AQI) rating (“Good”, “Moderate”, etc.) associated with 
PM 2.5 concentrations is shown with the vertical color bar. The pink and purple diamonds 
indicate fire locations as described for Figure 1. The 850 mb wind field vectors at 20Z 
from the Eta Data Assimilation System (ED AS) are plotted as arrows to show wind 
direction and speed. Winds at this level illustrate the flow in the lower troposphere. 
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Often this can be used to qualitatively show areas of convergence and divergence, and is 
one indicator of vertical air motion. An example wind vector is given in the legend. 
However, the magnitude and direction of this vector is only valid at the location given 
due to the map projection (Lambert’s conformal conic projection) used to display the 
data. This data fusion visualizes the relationship between the MODIS AOD and COT, 
hourly PM2.5 mass concentration and the air quality index, providing a pseudo-synoptic 
view of aerosol events across North America. 

3.2 Coincidence requirements 

For the correlation analysis presented in section 3.3, the data pairs of satellite AOD and 
ground based PM 2.5 must be co-located in space and time. For every ground station, the 
10x10 km 2 MODIS AOD observations that include the longitude and latitude of the site 
are accumulated. The hourly surface PM 2.5 data are then linearly interpolated to the time 
of each MODIS observation. Only surface observations within plus or minus one hour 
are considered for possible temporal coincidences. 

3.3 Time series analysis 

In figure 3, an example of the MODIS AOD and surface PM2.5 time series for the month 
of September 2003 is presented. The September 2003 time series plots for every ground 
station are provided in Appendix B. The station name, Metropolitan Statistical Area 
(MSA), and station ID are reported in the figure. If the site is in a rural area, the site name 
is listed as “Not in an MSA”. 

Figure 3 shows the surface PM2.5 data for both 1 -hour (solid line) and 24 -hour (dashed 
line) averages reported at an hourly frequency. The left vertical axis is mass 
concentration of PM2.5 and the right vertical axis is MODIS AOD. Coincident values are 
represented by symbols: 

• - MODIS AOD 

* - hourly PM2.5 mass concentration 

A- 24 hour average PM 2.5 mass concentration. 

Correlations are reported are for both 1-hour and 24 -hour average surface PM2.5 data. The 
correlations are derived from coincident MODIS AOD and PM2.5 data pairs as described 
above. N p0S s, is the total number of MODIS viewing opportunities over the site and 
Nmodis corresponds to the number of passes that optical depth values could be 
determined (cloud-free passes, etc. . .). The number of coincident data pairs used to 
determine the correlation is reported as Ncorr. Figure 3 provides critical information 
needed to determine if the MODIS AOD is indicative of PM2.5 concentrations at or near 
the surface. The vertical distance between coincident MODIS and PM2.5 points (as 
plotted) reflects the AOD/PM2.5 ratio. This ratio is another factor to be considered in 
determining if the aerosol is well mixed and at or near the surface or if the aerosol is 
aloft. If the distance is relatively small then the ratio is close to 1 . 6 / 100 ., and an inference 
can be made that the aerosols may be near the surface. 
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Figure 3. Time-series and correlations between MODIS AOD and hourly l-hour/24- 
hour average ground PM 2.5 mass concentrations for site 170434002 in Naperville, IL 
during September 2003. 

The vertical axes of the time-series plots are scaled as shown in the figure above for most 
of the sites listed in Appendix B. However, clipping of the hourly PM 2.5 mass 
concentration reading can occur when the reading exceeds the value of 100. Where 
appropriate, to show MODIS AOD values that exceed 1 .6 or 24-hour average PM 2.5 mass 
concentration values that exceed 100 , the vertical axes are proportionally expanded and 
labeled “***EXPANDED RANGE***”. 

3.4 Site-by-site correlation analysis 

Figure 4 summarizes the MODIS AOD and 1-hour PM 2.5 correlations derived from the 
September time series for each ground station across the United States and parts of 
Canada (see Appendix B). The size of the point plotted indicates the number of 
coincidences between MODIS AOD and hourly PM 2.5 concentrations during this period. 
The significance of the correlation generally increases with increasing number of 
coincidences. The color indicates the value of the correlation coefficient. This 
correlation summary provides a site specific and geographical perspective on how well 
the MODIS AOD observations depict the variability in surface PM 2.5 measurements. 
During this month, higher correlations are generally found in the eastern half of the US 
and in parts of the Pacific Northwest. 
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20030901-20030930 Correlation between AIRNOW 1 -hour PM2.5and MOD IS AOD 
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Figure 4 . National summary plot of correlations between MODIS AOD and hourly 
PM 2 .5. 
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Distribution of AOD/PM2.5 Correlations 



Site-by-Site Correlation 


Figure 5. Histogram of the site-by-site correlations between coincident MODIS AOD 
and AIRNow PM 2.5 measurements during September, 2003. 

Figure 5 shows the histogram of the site-by-site correlations during September, 2003. The 
histogram includes correlations for a total of 309 AIRNow sites. The distribution of 
correlations peaks at 0.8 with 51% of the AIRNow sites having correlations greater than 
0.5 (orange to red in Figure 4) and 11% of the AIRNow sites having negative correlations 
(blue in Figure 4) during September, 2003. 
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4.0 National Satellite and In-Situ Comparisons 

4.1 Maps of 40 km binned mean MODIS AOD statistics 

Figure 6 shows a map of the mean MODIS AOD for September, 2003. To construct this 
map, all of MODIS AOD granules (10 km x 10 km spatial resolution) obtained for 
September 2003 are mapped onto the Eta 40 km Lambert Conformal grids (185 x 129 
points). The mean, the standard deviation and the counts of MODIS AOD at each grid 
point are derived from the regridded MODIS AOD. The areas in black are where no 
MODIS AOD is retrieved over a 40 km x 40 km grid for the entire month due to either 
high surface reflectance or cloudiness. 


Figure 6. Map of the mean 40 km binned MODIS AOD for September, 2003. 

The September 2003 mean MODIS AOD is strongly influenced by the smoke plume 
from the fires in the Northwestern US. These fires led to very high aerosol loading along 
the Montana/North Dakota/Canada border during early September resulting in mean 
AOD values over 0.4 in this region. The local plume from the Northwestern US fires is 
evident along the Washington/Oregon border and extending into Eastern Washington 
State. Long-range transport of this smoke plume contributed to the high mean AOD in 
the Great Lakes region during September 2003. These regions influenced by wildfire- 
enhanced aerosol correspond to areas where high correlations exist between MODIS 
AOD and surface PM 2 .5 (Figure 4). Aerosol loading is generally high over Alabama and 
Mississippi. Agricultural burning in the lower Mississippi River valley and Alabama (see 
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Figure 1) contributed to this regional aerosol burning. High aerosol loading along the 
Gulf coast may be influenced by local sea breeze circulations transporting sea salt 
aerosols inland. Aerosol loading is generally low in the central plains region. High 
surface reflectance in the Southwestern US and Mexico results in false MODIS AOD 
retrievals over the desert Southwest. 


20030901-20030930 MODIS AOD Standard Deviation 
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Figure 7. Map of the standard deviation of the 40 km binned MODIS AOD for 
September, 2003. 

Figure 7 shows a map of the standard deviation of the 40 km binned MODIS AOD for 
September, 2003. The standard deviation is dominated by the influence of fires in the 
Northwestern US. The standard deviation of the MODIS AOD in regions influenced by 
smoke transport from these fires is as large as the mean AOD (Figure 6) due to the 
relatively short (~10 days) duration of this event. The standard deviation of the 40 km 
binned MODIS AOD is significantly lower than the mean 40 km binned AOD over the 
desert Southwest, indicating that high surface reflectivity introduces systematic, not 
random biases in the MODIS AOD retrievals over desert regions. Localized regions of 
high standard deviation in California, Nevada, Arizona, and Mexico reflect real variations 
in AOD associated with major urban areas. 
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20030901-20030330 MODIS AOD Counts 
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Figure 8. Map of the number of MODIS AOD retrievals within each 40 km bin for 
September, 2003. 

Figure 8 shows the number of MODIS AOD retrievals used to obtain the 40 km binned 
statistics. In the desert Southwest, low counts occur in regions with high mean AOD that 
are caused by high surface reflectivity (Western Nevada, Southern California, Salt Lake 
City, northeastern Arizona, southeastern New Mexico and west Texas). Counts are lower 
over the Gulf of Mexico and Atlantic Ocean due to sun glint interference in the AOD 
retrievals. Counts are lower over the Pacific due to missing MODIS AOD granules and 
persistent low-level cloudiness. 
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4.2 Site-by-site mean statistics 

The amount of information that MODIS AOD can contribute to characterization of the 
mean spatial distribution of aerosols at the AIRNow sites is quantified by comparing the 
site-by-site mean and standard deviations of the AIRNow PM 2.5 and MODIS AOD. 
Figures 9a and 9b show the site-by-site distribution of mean AIRNow PM 2.5 and mean 
MODIS AOD for September 2003. The September 2003 weighted national mean is 
1 1.75pg/m 3 and 0.28 for AIRNow PM 2.5 and MODIS AOD, respectively. Both mean 
distributions show elevated (above their respective weighted national mean) aerosols in 
the Central and South Eastern US, Los Angeles and Salt Lake City with generally low 
aerosols (below their respective weighted national mean) in the Northeastern US and 
Western Washington. MODIS shows relatively high mean AOD for the AIRNow sites in 
Eastern Washington while the mean surface PM 2.5 measurements remain low in spite of 
the large aerosol loading associated with the fires in the Northwest. This difference can 
be understood since aerosol loading in this region was likely transported above the 
boundary layer and therefore not sampled by the AIRNow surface sites. Similar 
influences of lofted aerosol loading are evident in Eastern Michigan and the Minneapolis- 
St. Paul metropolitan area where MODIS measured relatively high mean AOD but the 
mean AIRNow PM 2.5 measurements are generally low. 

Figures 10a and 10b show site-by-site distribution of the standard deviation of the 
coincident AIRNow PM 2.5 and MODIS AOD measurements for September 2003. The 
September 2003 weighted national average standard deviation is 6 . 6 pg/m 3 and 0.16 for 
AIRNow PM 2.5 and MODIS AOD, respectively. Both PM 2.5 and AOD show relatively 
high (above the weighted national average) standard deviations throughout the central 
and south central US, consistent with the standard deviation of the 40 km binned MODIS 
AOD (Figure 7). The PM 2.5 measurements show very high (>25pg/m 3 ) standard 
deviations for an individual AIRNow site in Los Angeles that is not present in the AOD 
measurements. The standard deviation at this highly urbanized site reflects very localized 
variability and is significantly larger than the mean PM 2.5 at this site during September 
2003. The MODIS AOD shows relatively high standard deviations at sites in Eastern 
Washington, Chicago, and Alabama. These high standard deviations most likely reflect 
variability associated with lofted aerosol loading since they are not as pronounced in the 
surface PM 2.5 standard deviations. 
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Figure 9. Site -by-site distribution of mean AIRNow PM 2.5 (9a) and MODIS AOD (9b) 
for September 2003. 
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Figure 10. Site-by-site distribution of AIRNow PM 2.5 (10a) and MODIS AOD (10b) 
standard deviations for September 2003. 
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4.3 Regional spatial statistics 

The comparison of maps of AIRNow PM2.5 and MODIS AOD site -by-site mean and 
standard deviation statistics shows that MODIS AOD can provide useful qualitative 
information about the spatial distribution of mean surface PM2.5 during September 2003. 
However, there is significant site-to-site variation in the agreement between the mean 
PM2.5 and AOD measurements. To quantify the spatial information content in the mean 
MODIS AOD we consider the correlations between site -by-site means and standard 
deviations within each of the EPA regions. Figure 1 1 shows a map of the EPA regions. 
Only continental US regions were considered in this analysis. 
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Figure 11. Map of EPA Regions 1-10. 


Figure 12 is an example of the spatial correlations between means and standard 
deviations of MODIS AOD and PM 2 .5 measurements for each AIRNow site within EPA 
regions 4 and 5 during September 2003. These regions have been combined so that the 
spatial information content of the MODIS AOD within the high aerosol loading over the 
Central and South-Central US can be quantified. (Appendix C summarizes the spatial 
correlations between AOD and PM 2 .5 means and standard deviations for each individual 
EPA region.) The spatial correlation between the mean AOD and PM 2.5 within EPA 
region 4 and 5 is low (-0.25). The spatial correlation between the AOD and PM 2.5 
standard deviations within EPA region 4 and 5 is somewhat higher (~0.4). In general, the 
spatial correlations within EPA regions are significantly lower than the site-by-site 
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temporal correlations shown in Figure 4. For example, the EPA region 4 and 5 weighted 
average of the site-by-site temporal correlations shown in Figure 4 is 0.59 during 
September 2003, which is over twice as high as the mean spatial correlation. 
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Figure 12. Spatial correlations between means and standard deviations of MODIS AOD 
and PM 2.5 measurements for each AIRNow site within EPA regions 4 and 5. 
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5.0 Conclusion 


This technical evaluation of the relationship between AIRNow surface PM 2.5 and 
satellite-observed AOD for a pilot study period of September 2003 indicates that MODIS 
AOD provides a daily, national perspective on atmospheric column aerosol loading 
which complements the AIRNow network by filling in gaps in the surface PM 2.5 
network. 

In addition, changes from day to day in the national perspective (MODIS AOD) provide 
insight into transport of aerosol, both laterally across the continent, and vertical ascent or 
descent. Good quantitative correlations with MODIS AOD and AIRNow surface PM 2.5 
are found when the aerosol is mostly in the boundary layer. Poorer quantitative results 
occur when the aerosol is lofted. Consequently, we have recommended that the MODIS 
Science Team devise a parameter to indicate the approximate vertical location of the 
aerosol (likelihood that the aerosol is near the surface or lofted) based on the meteorology 
and physics of the retrieval. 

Generally high agreement between MODIS AOD and AIRNow PM2.5 is found for the 
hourly data. However, MODIS AOD provides only qualitative information about the 
mean spatial distribution of surface PM2.5 during September 2003. The discrepancy 
between the generally high temporal correlations and generally low mean spatial 
correlations arises because the spatial statistics are influenced by site-to-site variations in 
the aerosol composition and the altitude of the aerosol loading, both of which lead to 
different mean AOD/ PM2.5 ratios. When these different mean PM2.5 /AOD ratios are 
combined into regional spatial correlations, the overall correlation decreases. In contrast, 
site-by-site background aerosol composition appears to be relatively constant during 
September 2003, so that significant events like the long-range transport of the smoke 
from the Northwestern wild fires result in episodic increases in both surface PM2.5 and 
AOD which are strongly correlated in time, even though individual sites have a different 
background relationship between AOD and PM2.5. 
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Appendix A 


Daily Satellite and EPA In-Situ Fusion National Maps 
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Appendix B 

Site-by-Site Satellite and EPA In-Situ Time Series 
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Appendix B contains Table Bl, Table B2, and the site-by-site satellite and in-situ time 
series plots for all of the ground stations used in this analysis. Table Bl lists pertinent 
information (EPA region, station ID, state, longitude, latitude, MSA number, MSA 
description, station name and monitor method) for the ground station sites in the United 
States organized by EPA region. In the case where the monitor method is left blank, that 
information was unavailable to the authors at the time of this printing. For the sites in 
Canada, only the station ID, province, longitude, latitude, and station name were 
available. The time series plots are listed in the same order as the Tables Bl and B2. 
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Table Bl. United States EPA Ground Station Sites 
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202090021 KS -94.64 39.12 3760 KANSAS CITY, MO-KS WY/KC FDMS-Grav 

295100085 MO -90.20 38.66 7040 ST. LOUIS, MO-IL BLAIR STREET TEOM Gravimetric 30 deg 

080010006 CO -104.94 39.83 2080 DENVER, CO Commerce City TEOM Gravimetric 30 deg 

080310002 CO -104.99 39.75 2080 DENVER, CO CAMP TEOM Gravimetric 30 deg 



060590001 CA -117.91 33.82 5945 ORANGE COUNTY, CA Anaheim 

060650012 CA -116.87 33.93 6780 RIVERSIDE-SAN BERNARDINO, Banning-South Hathaway Street 


B-12 



Region Station ID State Lon Lat MSA # MSA description Station Name Monitor Method 

9 061072002 CA -119.29 36.33 8780 VISALIA-TULARE, Visalia-N. Church Street 

PORTERVILLE, CA 
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Table B2. Canada Ground Station Sites 


Station ID 

Province 

Lon 

Lat 

Station Name 

000040103 

NB 



FREDERICTON 

000040203 

NB 

-66.01 

45.31 

FOREST HILLS 

000040901 

NB 

-67.08 

45.09 

ST ANDREWS 

000050105 

QC 

-73.58 

45.50 

Drummond 

000050110 

QC 

-73.64 

45.59 

Montreal-Nord 

000050126 

QC 

-73.93 

45.43 

STE-ANNE DE BELLEVU 

000050128 

QC 

-73.75 

45.47 

AEROPORT DE MONTREAL 

000050129 

QC 

-73.57 

45.66 

RIVI C8RE-DES-PRAIRIE 

000050131 

QC 

-73.75 

45.47 

Hochelaga 

000050308 

QC 

-71.22 

46.82 

DES SABLES 

000050801 

QC 

-72.54 

46.35 

Ursulines 

000051201 

QC 

-72.74 

46.55 

SHAWINIGAN 

000052101 

QC 

-73.64 

45.20 

L’ACADIE 

000052301 

QC 

-74.48 

46.03 

SAINT-FAUSTIN 

000054401 

QC 

-74.28 

45.12 

ST-ANICET 

000054501 

QC 

-73.44 

45.81 

L AS SOMPTION 

000054703 

QC 

-72.43 

46.35 

Becancour 

000060104 

ON 

-75.68 

45.43 

Ottawa 

000060204 

ON 

-83.04 

42.31 

Windsor Downtown 

000060302 

ON 

-76.51 

44.23 

Kingston 

000060424 

ON 

-79.39 

43.66 

Toronto Downtown 

000060512 

ON 

-79.86 

43.26 

Hamilton Downtown 

000060707 

ON 

-84.35 

46.53 

Sault Ste Marie 

000061004 

ON 

-82.41 

42.98 

Sarnia 

000061104 

ON 

-78.35 

44.30 

Peterborough 
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ON 

-74.74 

45.03 

Cornwall 
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ON 

-79.24 

43.16 

St. Catharines 
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ON 

-80.50 

43.44 

Kitchener 
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ON 

-79.45 

46.32 
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ON 

-82.22 

42.24 

Merlin 

000062501 

ON 
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44.30 
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ON 

-80.27 
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ON 

-78.93 
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ON 
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43.33 

Grand Bend 
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ON 
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Barrie 
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ON 

-80.04 

45.34 

Parry Sound 
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ON 

-81.16 

42.67 

Port Stanley 
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ON 

-77.40 

44.15 

Belleville 

000110119 

BC 

-122.98 

49.22 

Burnaby South 

000111003 

BC 

-122.23 

49.04 

Abbotsford 
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BC 
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49.16 

Chilliwack 
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BC 
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49.10 
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Appendix D 
Acronyms 


AOD - Aerosol Optical Depth 
AQI - Air Quality Index 
ARL - Air Resources Laboratory 
BAM - Beta Attenuation Monitor 

CIMSS - Cooperative Institute for Meteorological Satellite Studies 

COT - Cloud Optical Thickness 

ED AS - Eta Data Assimilation System 

EPA - United States Environmental Protection Agency 

ESA - Earth Science Applications 

ESE - Earth Science Enterprise 

FEM - Federal Eqivalent Monitor 

FRM - Federal Reference Monitors 

GOES - Geostationary Operational Enviro nm ental Satellite 
GRIB - Gridded Binary 

IDEA - Infusing satellite Data into Environ m ental Applications 

MODIS - Moderate Resolution Imaging Spectroradiometer 

MSA - Metropolitan Statistical Area 

NAMS - National Ambient Monitoring Stations 

NASA - National Aeronautics and Space Administration 

NCEP - National Centers for Environmental Prediction 

NESDIS - National Environmental Satellite, Data and Information Services 

NOAA - National Oceanic and Atmospheric Administration 

PM - Particulate Matter 

PM 2.5 - Particulate Matter less than 2.5 pm in diameter 
SLAMS - State and Local Air Monitoring Stations 
TEOM - Tapered Element Oscillating Microbalance 
US - United States 

WF ABBA - Wild Fire Automated Biomass Burning Algorithm 
WRI - World Resources Institute 


D-l 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and 
Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person 


shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 
01- 07-2004 

2. REPORT TYPE 3. DATES COVERED (From - To) 

Technical Memorandum 

4. TITLE AND SUBTITLE 

Comparison of Satellite Observations of Aerosol Optical Depth to Surface 
Monitor Fine Particle Concentration 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Kleb, Mary M.; Al-Saadi, Jassim A.; Neil, Doreen O.; Pierce, Robert B.; 
Pippin, Margaret R.; Roell, Marilee M.; Kittaka, Chieko; and Szykman, 
James J. 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

23-613-23-01 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-2199 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L-19039 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSOR/MONITOR'S ACRONYM(S) 

NASA 

11. SPONSOR/MONITOR'S REPORT 
NUMBER(S) 

NASA/TM-2004-2 13248 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 45 

Availability: NASA CASI (301) 621-0390 Distribution: Standard 

13. SUPPLEMENTARY NOTES 

An electronic version can be found at http://techreports.larc.nasa.gov/ltrs/ or http://ntrs.nasa.gov 


14. ABSTRACT 

Under NASA’s Earth Science Applications Program, the Infusing satellite Data into Environmental Applications (IDEA) project examined the relationship 
between satellite observations and surface monitors of air pollutants to facilitate a more capable and integrated observing network. This report provides a 
comparison of satellite aerosol optical depth to surface monitor fine particle concentration observations for the month of September 2003 at more than 300 
individual locations in the continental US. During September 2003, IDEA provided prototype, near real-time data-fusion products to the Environmental 
Protection Agency (EPA) directed toward improving the accuracy of EPA’s next-day Air Quality Index (AQI) forecasts. Researchers from NASA Langley 
Research Center and EPA used data from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument combined with EPA ground network data 
to create a NASA-data-enhanced Forecast Tool. Air quality forecasters used this tool to prepare their forecasts of particle pollution, or particulate matter less 
than 2.5 microns in diameter (PM2.5), for the next-day AQI. The archived data provide a rich resource for further studies and analysis. The IDEA project uses 
data sets and models developed for tropospheric chemistry research to assist federal, state, and local agencies in making decisions concerning air quality 
management to protect public health. 


15. SUBJECT TERMS 

Particulate matter; Air quality forecast; MODIS; AIRNow; PM2.5; Aerosol; Air pollution; Earth Science data applications; Satellite data 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 

U 

u 

u 

uu 


18. NUMBER 
OF 

PAGES 


19a. NAME OF RESPONSIBLE PERSON 

STI Help Desk (email: help@sti.nasa.gov) 


19b. TELEPHONE NUMBER (Include area code) 

(301) 621-0390 


Standard Form 298 (Rev. 8-98) 

Prescribed by ANSI Std. Z39.18 


278 



































